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We review recent experimental and theoretical developments in inclu- 
sive semileptonic B — )■ Xciiy decays. In particular, we discuss the deter- 
mination of \Vcb\ and of the heavy quark masses through fits based on the 
Operator Product Expansion. 
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Semileptonic b ^ c transitions play an important role in various aspects of heavy 
flavor physics. They allow for a determination of the CKM matrix element Vet, which 
is a crucial input in all Unitarity Triangle analyses: for instance the (p, fj) constraint 
determined by Ek is very sensitive to the precise value of the Wolfenstein parameter 
A, which is essentially determined by \Vcb\- Moreover, the bottom quark mass and 
the hadronic parameters extracted from fits to inclusive semileptonic and radiative 
moments are key inputs for the inclusive \Vub\ determination, the normalization of 
rare B decays like B — )■ X^'j, and various other B physics applications. 

We will review the present theoretical and experimental status of inclusive B 
Xc^v decays, with particular emphasis on recent developments, among which new 
calculations and measurements. 

1 Theoretical framework 

Our understanding of inclusive semileptonic B decays rests on a simple idea: as all 
final states are summed over in inclusive decays, the final quark hadronizes with unit 
probability and the transition amplitude is sensitive only to the long-distance dy- 
namics of the initial B meson, which can indeed be factorized. An Operator Product 
Expansion (OPE) allows us to express the non-perturbative physics in terms of ma- 
trix elements of local operators of dimension (i > 5, while the Wilson coefficients are 
perturbative [U |2]. The leading term in this double expansion in as and Aqco/'^fe is 
given by the free h quark decay, and the first corrections are 0{as) and 0(AQCj~|/m^). 
The relevant parameters are the heavy quark masses mb and mc, the strong coupling 
as, and the matrix elements of the local operators: p^ and p^. at 0(l/m^), pf, and p\g 
at 0(l/mf), etc. Since the OPE is valid only for sufficiently inclusive measurements 
and away from perturbative singularities, the relevant quantities to be measured are 
global shape parameters (the first few moments of various kinematic distributions) 
and the total rate. The former give information on the masses and matrix elements, 
the latter on \ Vcb\- The OPE parameters describe universal properties of the B meson 
and of the quarks and are useful in many applications. 

The main ingredients for an accurate analysis of the experimental data have been 
known for some time. Two implementations are currently employed, based on either 
the kinetic scheme [H HI [5] or the IS* scheme [6]. They both include terms through 
0(a^/3o) [Z] and 0{l/m\) |S] but they use different perturbative schemes, include 
a somewhat different choice of experimental data under specific assumptions, and 
estimate the theoretical uncertainty in two distinct ways. Nevertheless, it is reassuring 
that, as we will show below, the two methods yield very close results for \Vcb\- 

The complete two-loop perturbative corrections to the width and moments of the 
lepton energy and hadronic mass distributions have been recently computed [H HU] by 
both numerical and analytic methods. In general, using as{frtb) in the on-shell scheme. 



1 



the non-BLM corrections amount to about —20% of the two-loop BLM corrections. 
In the kinetic scheme with cutoff /x = IGeV, the perturbative expansion of the total 
width is 

r[B^ X,eu] oc 1 - 0.96 - 0.48/3o f-V + 0.82 f — V + O(a^) ^ 0.916 (1) 

Higher order BLM corrections to the width and moments are also known [3l [7j. The 
resummed BLM result is numerically very close to the NNLO one [3]. The residual 
perturbative error in the total width is therefore about 1%. 

Since the numerical results of are available for a variety of lepton energy 
cuts and values of rric/mb, it is now possible to implement them in a global fit. In 
the normalized leptonic moments the perturbative corrections cancel to large extent, 
independently of the scheme, as hard gluon emission is comparatively suppressed. 
This pattern of cancellations, crucial for a correct estimate of the theoretical error, 
is confirmed by the complete 0{a1) calculation, although the numerical precision of 
the available results is not always sufficient to improve the final accuracy. The actual 
implementation in the kinetic scheme is under way. 

Another source of significant theoretical uncertainty are the 0{asAQQj^/ml) cor- 
rections to the width and to the moments. Only the ©(a^/i^/m^) terms are known [TT] . 
A complete calculation of these effects has been recently performed in the case of in- 
clusive radiative decays [T2j, where the 0{as) correction increase the coefficient of fiQ 
in the rate by almost 20%. The extension of this calculation to the semileptonic case 
is in progress. In view of the numerical importance of 0(l/mf) corrections, if the 1% 
precision in the width is to be reached, the effects 0{as/ml) will also be necessary. 

For what concerns the higher order power corrections, a thorough analysis of 
0{l/mf) and ©(l/mg) effects has just appeared [13]. The main problem is the pro- 
liferation of non-perturbative parameters: e.g. as many as nine new expectation values 
appear at 0{l/mf). As they cannot be fitted from experiment, in Ref. [13] they are 
estimated in the ground state saturation approximation, reducing them to the known 
0{l/ml'^) parameters. In this approximation the total 0(1 /mg^) correction to the 
width is about +1.3%. The O^I/itiq) effects are dominated by 0{l/mlml) Intrinsic 
Charm contributions, amounting to -|-0.7% [H]. The actual effect on |Vcfe| depends 
also on the corrections to the moments. The authors of [13] estimate that the overall 
effect on \Vcb\ is a 0.4% increase, consistent with our preliminary implementation of 
these effects in the kinetic global fit. While this sets the scale of higher order power 
corrections, it is yet unclear how much the result depends on the assumptions made 
on the expectation values. 

The first two moments of the photon energy distribution in B Xg'y are also rou- 
tinely included in the semileptonic fit. They are sensitive to and /i^ in particular. 
However, the experimental lower cut on the photon energy introduces a sensitivity to 
the Fermi motion of the 6-quark inside the B meson and tend to disrupt the OPE. 
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One can still resum the higher-order terms into a non-local distribution function and 
since the lowest integer moments of this function are given in terms of the local OPE 
parameters, one can parameterize it assuming different functional forms [3] . Another 
serious problem is that only the leading operator contributing to inclusive radiative 
decays admits an OPE. Therefore in principle unknown 0{asA/mb) contributions 
should be expected [15] and radiative moments should be considered with care in the 
context of high precision analyses. 

2 Measurements of moments 

BaBar has recently published a study of the hadronic mass spectrum mx in inclusive 
decays B — > Xcii' [IS]. The main steps of this analysis, based on a data sample of 
232 million T{4S) — )■ BB events, are: First, the decay of one B meson in the event is 
fully reconstructed in a hadronic mode (-Stag) and the associated tracks and clusters 
are removed from the event. Such a i?tag candidate can be found in about 0.4% of 
the T(4S') events with a signal purity of about 80%. Then, the semileptonic decay 
of the second B meson in the event (-Bsig) is selected by searching for an identified 
charged lepton (electron or muon) with momentum above 0.8 GeV/c. Finally, all 
remaining particles in the event are combined to reconstruct the hadronic X system. 
The resolution in mx is improved by a kinematic fit taking into account 4-momentum 
conservation and the consistency of the missing mass with a zero mass neutrino. 

Still, the observed mx spectrum is distorted by resolution and acceptance effects 
and cannot be used directly to obtain the hadronic mass moments. BaBar implements 
a linear correction to obtain the true moments from the reconstructed ones. Different 
corrections are applied depending on the X system multiplicity, -Emiss — cpmiss and the 
lepton momentum. In this way, BaBar measures the moments of the hadronic mass 
spectrum up to (m,^) for minimum lepton energies ranging between 0.8 and 1.9 GeV. 

This study also updates the previous BaBar measurement of the lepton energy mo- 
ments in 5 — !■ Xcii' [13 using new branching fraction measurements for background 
decays and improving the evaluation of systematic uncertainties. Also, the first mea- 
surement of combined hadronic mass and energy moments (n^) with k = 2,4,6 is 
presented, where the latter are defined as n\ = m^c^ — 2AEx + A^, with mx and 
Ex the mass and the energy of the X system and A a constant fixed to 0.65 GeV. 

BaBar interprets their data using the OPE in the kinetic scheme [21 IH IS] and 
performs a simultaneous fit to 12 hadronic mass moments (or 12 combined mass- 
energy moments), 13 lepton energy moments (including partial branching fractions 
as 'zero order' moments), and 3 photon energy moments in 5 — )■ Xg'j [HI [19]. The 
results are given in Table [TJ 

Also the Belle collaboration has obtained measurements of the lepton energy Ei 
and the hadronic mass spectrum mx in B ^ Xc^v using 152 million T(4S') — )■ 
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Table 1: Results of the OPE fits in the kinetic scheme to the BaBar data \T6\. The 
first uncertainty quoted is experimental, the second theoretical. 





Hadronic moments 


Mass-energy moments 


\Vcb\ (10-3) 


42.05 ± 0.45 ± 0.70 


41.91 ±0.48 ±0.70 


rub (GeV) 


4.549 ±0.031 ±0.038 


4.556 ±0.034 ±0.041 


B{B XJu) (%) 


10.64 ±0.17 ±0.06 


10.64 ±0.17 ±0.06 


xVndf. 


10.9/28 


8.2/28 



Table 2: Results of the OPE fits in the kinetic and IS schemes to the Belle data [23] 





Kinetic scheme 


IS scheme 


\Vcb\ (10-3) 


41.58 ±0.90 


41.56 ±0.68 


B{B XJu) (%) 


10.49 ±0.23 


10.60 ±0.28 


xVndf. 


4.7/18 


7.3/18 



BB events [201 EI]- The experimental method is similar to the BaBar analysis dis- 
cussed previously, i.e., one B meson is fully reconstructed in a hadronic mode and 
an identified lepton is required to select semileptonic decays of the second B. In 
the Belle analyses acceptance and finite resolution effects in the Ei and mx spectra 
are corrected by unfolding using the SVD algorithm [22j. Belle measures (Ef) for 
k = 0,1, 2, 3, 4 and minimum lepton energies ranging between 0.4 and 2.0 GeV. Mo- 
ments of the hadronic mass (rrix) sue measured for = 2,4 and minimum lepton 
energies between 0.7 and 1.9 GeV. 

To obtain \Vcb\, Belle fits 14 moments of the lepton energy spectrum, 7 hadronic 
mass moments and 4 moments of the photon energy spectrum in i? — )■ Xs'j [23] to 
OPE expressions derived in the kinetic [3l HI [5] and IS schemes [6|. Both theoretical 
frameworks are considered independently and yield very consistent results with the 
Belle data. Table |2} 

3 Global HFAG fit 

The Heavy Flavor Averaging Group (HFAG) has performed as global analysis of in- 
clusive observables in _B — t- Xc£z/ and B — )■ Xs'j decays to determine the 6-quark 
mass rrib and the higher order parameters in the OPE description of these decays. 
This analysis combines data from the BaBar, Belle, CLEO, CDF and DELPHI ex- 
periments. 
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Table 3: Experimental data used in the HFAG analysis of inclusive B — )■ Xcii' and 
B — )• decays. In the table, (-E^), (mx) and (Ej^) refer to the moments of the 
lepton energy and hadronic mass spectrum in i? — )■ Xcii^ and to the photon energy 
moments in 5 — )■ Xgj, respectively. The index k specifies the order of the moments 
used. 



BaBar 


{E^): fc = 0,l,2,3 [ISIIIZI, (m|): A; = 2, 4, 6 [IS], (El^): k = 1,2 


Belle 


{E^y. k = 0,1,2,3 mAm'k)- A; = 2,4 [21], (E^^): A; = 1,2 [25] 


CDF 


(m^): A; = 2,4 pS] 


CLEO 


(m^): A; = 2,4 [27], (E!;): k = l^ 


DELPHI 


{E^): A; = 1,2,3 [29], (m^): A; = 2,4 [29] 



Table 4: Results of the HFAG global fit in the kinetic scheme. The first error on 
\Vcb\ is the uncertainty from the global fit, the second is the error in the average 
B lifetime and the third error is an additional theoretical uncertainty arising from 
the calculation of \Vcb\- 



Input 






(10-3) m^i^ (GeV) fil (GeV^) xVndf. 


all moments 
Xciiy only 


41.85 ±0.42 ±0.09 ±0.59 4.591 ± 0.031 0.454 ± 0.038 29.7/59 
41.68 ±0.44 ±0.09 ±0.58 4.646 ± 0.047 0.439 ± 0.042 24.2/48 



The global fit is done both with expressions derived in the kinetic [3l HI |5] and 
IS schemes In both cases 7 free parameters are determined in the fit. The only 
external input used in the analysis is the average B^ and B^ lifetime. 

The data used in the global fit is listed in Table |3j In total 66 measurements - 
29 from BaBar, 25 from Belle and 12 from other experiments - are used. Note that 
the analysis in the IS scheme still uses the BaBar 2004 hadronic moment measure- 
ments |24] . 

The results of the global fit in the kinetic scheme are given in Table |4] and Fig. [Tj 
The results of the IS scheme analysis are shown in Table |5j In both cases, the results 
with all moments and with B — t- Xciiy moments only are quoted. 

4 Discussion 

We have seen that the fits discussed in the previous section determine nib quite 
precisely. How does this rrib determination compare with alternative determina- 
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m, (GeV) m, (GeV) 

Figure 1: A^^ = 1 contours for the HFAG global fit in the kinetic mass scheme. 



Table 5: Results of the HFAG global fit in the IS scheme. 



Input 






(10-^) ml^ (GeV) Ai (GeV^) x'/ndf. 


all moments 
Xcd-v only 


41.87 ±0.25 4.685 ±0.029 -0.373 ± 0.052 32.0/57 
42.31 ±0.36 4.619 ±0.047 -0.427 ± 0.057 24.2/46 



tions [301 ED [32]? Semileptonic moments do not measure rri}, well. As illustrated 
in Fig. [2| they rather identify a strip in the (rric, mi,) plane along which the minimum 
is quite shallow, and |V"cb| basically constant (straight lines). The global kinetic fit 
selects an (mc, rrih) region compatible with the loose PDG-2007 [33] bound^ as well 
as with the precise e^e~ sum-rules determinations [301 |3T], of course after conversion 
to the kinetic scheme. This conversion is known to 0{a1) and entails a non-negligible 
error, of about 40 MeV for the conversion from m^^(mfe) to m^*"(lGeV) and about 
10 MeV for that from mp(3GeV) to m^^"(lGeV). 

It turns out that the semileptonic fit, and in particular its determination of the 
masses and the other OPE parameters, is very sensitive to various details. For in- 
stance, the assumptions on the correlations between theoretical errors for moments 
evaluated at different cuts have a clear impact on the OPE parameters, while the 
value of \Vcb\ remains quite stable. Such theoretical correlations are obviously hard 
to estimate. The present HFAG fit follows the procedure outlined in [M], assuming 
100% correlation between moments calculated at different values of Ecuti the lower 

^ Later editions of PDG have stretched the uncertainties in an abnormal way. 
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Figure 2: Different charm and bottom quark: determinations in the kinetic mass 
scheme. The elhpses represent the PDG-2007 ranges (large green), a global semilep- 
tonic fit that differs slightly (see text) from the HFAG one (red), the Karlsruhe (pink) 
and Hoang et al. (blue) sum-rules determinations. 

cut on the lepton energy. This very strong assumption distorts the fit, leading to high 
values of m^^c, even outside the PDG range, and to underestimating the uncertainty 
of all non-perturbative parameters. On the opposite extreme, no correlation between 
close values of Ecut is unreasonable. A more realistic approach, adopted in the fit 
shown in Fig. 2, consists in taking into account the Ecut dependence and correlations 
of the known OPE calculation. It leads to slightly lower rrib^c with larger errors. A 
detailed discussion will be presented elsewhere [35] . 

A related question concerns the role of radiative moments in the fits: as shown 
above they help fixing mi,. But the fit is almost identical if one replaces them with the 
loose bound m^^(mfe) = 4.20(7) GeV given by PDG in 2007. Indeed, the inclusion 
of external, well-founded constraints in the fit can be very useful: it decreases the 
errors and neutralizes the potential weight of theoretical correlations. As semilep- 
tonic decays do determine precisely a linear combination of mft^c, a way to maximally 
exploit their potential consists in fitting directly m^^(3GeV) instead of the kinetic 
charm mass (this is possible and avoids the scheme conversion error), and including 
in the fit one of the recent very precise nic determinations. As an illustration we 
have used m^^(3GeV) = 0.986(13) GeV by the Karlsruhe group [30], and obtained 

"(IGeV) = 4.535(21) GeV, which translates into mf^{mb) = 4.165(45) GeV. This 
value for the bottom mass is perfectly consistent with the Karlsruhe group's own rrib 
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determination, m^^(mfe) = 4.163(16) GeV. The results of Refs. [3T] and |32] are also 
consistent. 

The kinetic scheme fitting routines are now undergoing a major upgrade, concern- 
ing the inclusion of higher order effects, the possibility to change the perturbative 
scheme, and the inclusion of additional constraints in the fit. The preliminary results 
we have just shown indicate that an uncertainty of about 20 MeV on rub can be reli- 
ably reached if an independent, precise determination of rric is employed. In view of 
this progress and of the calculations recently completed or under way, we believe that 
a 1% determination of \Vcb\ can be reached, although some work is still necessary. 
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